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We studied the spin torque switching of a single free layer in the thermally activated region by numerically solving the 
Landau-Lifshitz-Gilbert equation. We found that the temperature dependence of the switching time of the in-plane mag- 
netized system is nonlinear, which means b ± 1. Here, b is the exponent of the current term in the switching rate formula 
and has been widely assumed to be unity. This result enables us to evaluate the thermal stability of spintronics devices. 



Thermal stability is an important property of ferromagnetic 
materials for spintronics device applications such as spin ran- 
dom access memory (Spin RAM) and microwave oscillator. 
For example, the greater the thermal stability, the longer the 
Spin RAM retention time. 12) Experimentally, the thermal sta- 
bility is evaluated by measuring the magnetization switching 
of a free layer in the thermally activated region, and analyz- 
ing the time evolution of the switching probability with the 
formula, 3 ' 



- f v(f')df' 
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P(t) = 1 - exp 
where v(f) is the switching rate given by 

#appl 



v(f) = / exp 



Ao 1 + 



1 - 



(1) 



(2) 



Here, / is the attempt frequency. The thermal stability, Ao, 
of a uniaxially anisotropic ferromagnet is defined as Aq = 
MH K V/(2k B T), where M, H K , V, and T are the magnetiza- 
tion, uniaxial anisotropy field, free layer volume, and temper- 
ature, respectively. 7 C is the critical current of the spin torque 
switching at zero temperature. The exponent of the current 
term is denoted as b. Equations (1) and (2) were analytically 
derived for the uniaxially anisotropic system by solving the 
Fokker-Planck equation, and b — 2 as shown in refs. 4-6 ' On 
the other hand, for an in-plane magnetized system, which has 
easy and hard axes along and normal to the film plane, re- 
spectively, and does not have axial symmetry, it is difficult to 
derive the analytical formula of the switching rate from the 
Fokker-Planck equation. However, since eq. (2) is the general 
form of the switching rate following the Arrhenius law, 7 ' eq. 
(2) has been widely used to determine the thermal stability of 
the in-plane magnetized system. 1 ' 2 ' 

Let us discuss the differences between the uniaxially 
anisotropic and in-plane magnetized systems from the stand- 
point of the Fokker-Planck approach. In the uniaxially 
anisotropic system, the magnetization dynamics shown in Fig. 
1(a) is described by one variable (the angle from the easy 
axis). Moreover, the effect of the spin torque can be included 
in the effective potential, 4 ~ 6 ' and can be regarded as an addi- 
tional term to the applied field. Then, as the field switching 
problem, 3 ' the switching rate formula can be obtained analyt- 
ically. On the other hand, in the in-plane magnetized system, 
the magnetization dynamics shown in Fig. 1(b) is described 
by two angles (the angles from the easy and hard axes). Also, 
no effective potential can be introduced to describe the spin 
torque effect. These make it difficult to calculate the switch- 
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Fig. 1 . (a) Schematic views of the magnetization switching in the uniaxi- 
ally anisotropic system. The black arrow and red line represent the magneti- 
zation and its orbit, respectively, (b) Magnetization switching in the in-plane 
magnetized system. 

ing rate analytically. 

In 2004, based on the Landau-Lifshitz-Gilbert (LLG) equa- 
tion, Koch et aZ. 8 ' showed that b — 1 for the in-plane mag- 
netized system. After that, eq. (1) with b — 1 has been 
widely used to analyze experiments. 1 ' 2 ' Recently, however, 
we 9 ' pointed out that the result in ref. 8 ' should be regarded 
as an approximate one of the uniaxially anisotropic system 
due to the assumption used in their calculation, and is not ap- 
plicable to the in-plane magnetized system. Then, the natu- 
ral question becomes: What is the value of b for the in-plane 
magnetized system? It should be noted that the value of b 
significantly affects the evaluation of the quality of the spin- 
tronics devices. 9 ' For example, the retention time of the Spin 
RAM estimated by using the theory with b = 2 is several or- 
ders of magnitude longer than that estimated by using b = 1 . 
Thus, the determination of the value of b is important for spin- 
tronics. In particular, since the in-plane magnetized system is 
conventionally used as the fundamental structure of the Spin 
RAM, the determination of b in this system is an attractive 
problem for practical application. 

In this letter, we studied the spin torque switching of the 
single free layer in the thermally activated region by numeri- 
cally solving the LLG equation. According to our recent stud- 
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10, 11) 



investigation of the temperature dependence of the 
switching time enables us to estimate the value of b. By com- 
paring these previous works with the current results, we found 
that b ~ 3 for in-plane magnetized system. This value may 
not be a universal one. However, the important point is that 
our results indicate b + 1 . 

Before proceeding to the in-plane magnetized system, we 
first study the switching of the uniaxially anisotropic system 
to show the consistency between the LLG and Fokker-Planck 
approaches for the spin torque system. As mentioned above, 
for the Fokker-Planck approach to the uniaxially anisotropic 
system, the spin torque can be regarded as an additional 
term to the applied field. Thus, the switching probability for 
= {h, 0) should be identical to that for {\h\, = (0, h), 
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Fig. 2. Time evolutions of (a) the averaged m z and (b) the switching prob- 
ability under the effect of the applied field only. The system is uniaxially 
anisotropic. The time evolutions of (c) the averaged m z and (d) the switching 
probability under the effect of the spin torque only are also shown. 

where < h < 1, h = H app \IH^, and j = -I/I c : see eq. (2). 
We used this fact to check the consistency. 

The details of the calculation are as follows. We assumed 
that the magnetization dynamics is described by the LLG 

equation: 3,12 ' 13 - 1 

dm dm 

— = -ym x H + yH s m x (p x m) -ymxh + o-mx — , (3) 
at at 

where m, H, y, and a are the unit vector pointing in the direc- 
tion of the magnetization, the magnetic field, the gyromag- 
netic ratio, and the Gilbert damping constant, respectively. 
H s = HT]I/{2eMV) and p are the strength of the spin torque 
and the unit vector pointing in the direction of the magneti- 
zation of the pinned layer, respectively. The positive current 
with the spin polarization r\ is defined as the electron flow 
from the pinned to the free layer. The random field satisfies 
{hi{i)hj(t')) = [2ak B T KjMV)]S u S{t-t'), where/,;' = x,y,z. 3) 
In eq. (3), we use the macrospin model because the current- 
resistance curve in ref. '' shows rapid changes of the resistance 
between the parallel (m = p) and antiparallel (m = -p) align- 
ments, which means a uniform rotation of the magnetization. 
We assume that both the easy axis of the free layer and p are 
parallel to the z-axis (p = (0, 0, 1)). The initial state is set to be 
m(0) = (0, 0, 1). The 4 th -order Runge-Kutta method was em- 
ployed to solve eq. (3). The magnetization dynamics is aver- 
aged over 10 5 samples. The switching probability is obtained 
by counting the sample numbers in which m z < -0.9. Once 
the magnetization reaches m z < -0.9, we regard the sample 
as the switched system, even if the magnetization returns to 
the region m z > -0.9 due to the thermal fluctuations. 

For the uniaxially anisotropic system, the magnetic field 
is given by H = (0, 0, H app \ + H^m z ). The critical cur- 
rent from the parallel to the anti-parallel alignment is I c = 
-[2aeMV7(%)](i/ app i + H K ). According to the definition, the 
negative h and j favor the antiparallel alignment. Figures 2(a) 
and 2(b) show the time evolutions of the averaged m z and 
switching probability under the effect of the applied field 
only = 0), respectively. The values of the parameters are 
taken to be M = 1000 emu/c.c, H K = 200 Oe, y = 17.64 
MHz/Oe, a = 0.01, and T = 300 K. The thickness and cross- 



sectional area of the free layer are taken to be 2.5 nm and 
n x 80 x 35 nm 2 , respectively. 2 ' The magnitude of the ap- 
plied field, H. dpp \ - -180 Oe, is 90% of the anisotropy field 
(h = -0.9). On the other hand, in Figs. 2(c) and 2(d), the time 
evolutions of the averaged m z and switching probability under 
the effect of the spin torque only (h = 0) are shown, where 
j = -0.9. The switching probability for (\h\, \j\) = (0.9,0) is 
almost identical to that for (\h\, \ j\) = (0,0.9). This result in- 
dicates the consistency between the LLG and Fokker-Planck 
approaches for the spin torque system, and can be regarded as 
evidence of b — 2 from the LLG approach. 

Figures 2(b) and 2(d) show the waiting time during which 
the switching probability remains zero, while eq. (1) predicts 
a slight increase of the switching probability from t — 0. The 
waiting time was found in both experiments 14 ' and numeri- 
cal calculation of the Fokker-Planck equation. 5 ' Although the 
waiting time is important for device application, further study 
is beyond the scope of this letter. 

Next, we move to the in-plane magnetized system. The 
magnetic field is given by H = (0, -4nMm y , H- dpp \ + H-gjn^, 
where the y-component includes the demagnetization field. In 
this system, the value of b can be estimated as follows. By 
assuming the constant attempt frequency /o and the sweep 
current 7(f) = xt with the sweep rate x as is done in the exper- 
iments, 15 ' eqs. (1) and (2) with H dpp \ = reduce to 



P(t) = 1 - exp 
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where y(J3, z) is the lower incomplete F function. As shown in 
ref., 9 ' the probability density, dP{t)/dt, has its maximum at a 
certain time t . We call t the switching time, which is given by 
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for b = 1 , and 
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for b > 1 . Here, plog(z) is the product logarithm. By assuming 
that Ao oc l/T, in the low-temperature region, the temperature 
dependence of t is approximately linear if b — 1, whereas it is 
nonlinear if b > l. 1011 ' Thus, we investigated t by solving eq. 
(3), and compared it with eq. (5) or (6). 

In the calculation, the current sweep rate is taken to be 30 
A/s. The switching time is obtained by fitting the probabil- 
ity density with a Gaussian curve. It should be noted that the 
magnitude of the critical current, 7 C , in eq. (5) or (6) is usually 
larger than that estimated on the basis of the instability of the 
initial state, 16 ' -[2aeMV/(hrj)](H K + 2nM) = -0.54 mA {q is 
taken to be 0.8), because the instability of the initial state does 
not guarantee switching. Thus, the switching time shown be- 
low is longer than [2aeMV/(xhi])](H K + 2nM) = 18 fxs, and 
eq. (5) or (6) has four fitting parameters, I c , b, A^T, and fg. 

Figure 3 shows the time evolutions of the switching proba- 
bility and probability density at T = 1 K. The probability den- 
sity can be well reproduced by the Gaussian curve as indicated 
by the dotted line. The probability densities for various tem- 
peratures (< 20 K) and their switching time are shown in Figs. 
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Fig. 3. Time evolutions of the switching probability of the in-plane mag- 
netized system (black) and the probability density (red) at T = IK. The 
dotted line (blue) is the Gaussian fit, where the central time is defined as the 
switching time. 
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Fig. 4. (a) Probability densities for the various temperatures, (b) Temper- 
ature dependence of the switching time of the in-plane magnetized system. 
The inset shows the switching time in the very low temperature region. 

4(a) and 4(b), respectively. The switching time shows nonlin- 
ear temperature dependence, which means b + 1 . While the 
estimated values of AoT and /o depend heavily on the initial 
values of the fitting, the estimated value of b depends weakly 
on the initial value, and is b ~ 3. 

In an infinite demagnetization field limit, the switching is 
completely limited in the in-plane, i.e., one-dimensional mo- 
tion as the uniaxially anisotropic system, and b = 2. 6) In a 
realistic system, while there are infinite numbers of possi- 
ble switching paths, the averaged magnetization dynamics is 
determined mainly by the path along the in-plane: see Fig. 
1(b). However, because the demagnetization field is finite, the 
magnetization can move to the hard axis direction. Due to the 
contribution of such paths, b should deviate from 2. In other 
words, our result may indicate that b of the in-plane mag- 
netized system depends on the material parameters, such as 
the demagnetization coefficient. Although the investigation on 
such dependence is under consideration, the important point 



is that our results indicate b + 1 for the in-plane magne- 
tized system. We also briefly consider why b estimated above 
is larger than 2. Compared with the infinite demagnetization 
field limit, the realistic system has a large number of switch- 
ing paths, and thus, a high switching probability. Since the 
large b corresponds to a large number of switching events, b 
may be larger than 2. However, more accurate discussion re- 
quires the analytical study based on the Fokker-Planck equa- 
tion, and is beyond the scope of this letter. 

In summary, we studied the spin torque dependence of the 
magnetization switching probability by numerically solving 
the Landau-Lifshitz-Gilbert equation. We found that the tem- 
perature dependence of the switching time is nonlinear, i.e., 
b + 1, where b is the exponent of the current term in the 
switching rate formula. The results will be important for eval- 
uating the thermal stability of Spin RAM. 

The authors would like to acknowledge S. Yuasa, H. Kub- 
ota, Y. Suzuki, H. Sukegawa, and S. Mitani for the valuable 
discussions they had with us. 
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